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Abstract 

An investigation of the spin excitation spectrum of charge ordered (CO) a'- 
NaV205 is presented. We discuss several different exchange models which 
may be relevant for this compound, namely in- line and zig-zag chain models 
with weak as well as strong inter- chain coupling and also a ladder model 
and a CO/MV (mixed valent) model. We put special emphasis on the im- 
portance of large additional exchange across the diagonals of V- ladders and 
the presence of exchange anisotropics on the excitation spectrum. It is shown 
that the observed splitting of transverse dispersion branches may both be in- 
terpreted as anisotropy effect as well as acoustic- optic mode splitting in the 
weakly coupled chain models. In addition we calculate the field dependence 
of excitation modes in these models. Furthermore we show that for strong 
inter- chain coupling, as suggested by recent LDA-I-U results, an additional 
high energy optical excitation appears and the spin gap is determined by 
anisotropies. The most promising CO/MV model predicts a spin wave dis- 
persion perpendicular to the chains which agrees very well with recent results 
obtained by inelastic neutron scattering. 
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I. INTRODUCTION 



Transition metal- oxygen pyramids are ideal building blocks to obtain insulators with 
low D structures of 3d- ions like chains or ladders. Their locahzed spins exhibit collective 
quantum properties at low temperatures, e.g. spin gap formation in S—1 chains or S—1/2 
ladders. In addition there is the possibility of spin gap appearance due to the spin- Pcierls 
(SP) mechanism which causes dimerization of the chain. The standard example is CuGcOs^. 
Recently a'- NaV205 which has the Trellis lattice structure with alternating shifted ladders 
(Figs.l and 6) was investigated for similar reasons: observation of a superstructure below 
Tc— 33 K and subsequent spin gap formation as witnessed by a drop of the susceptibility 
below Tc^. However it is clear now that this compound does not exhibit a standard SP- 
transition because above Tc it is a homogeneous mixed valent (MV) insulator with one 3d- 
electron per V-V rung. Therefore above Tc a'- NaV205 is a quarter filled ladder system 
with equivalent V- sites instead of a family of half filled (atomic spin) chains. This was 
concluded from x- ray^~^ and NMR- experiments^^. They also show that below Tc in the 
dimerized state two inequivalent V- sites exist. Therefore a charge ordering (CO) transition 
which localizes the V 3d- electrons on one site of each rung of the ladders must take place 
at Tc. Possible CO- structures have been discussed by various authors*^'®' but so far the 
real low temperature structure remains controversial. In general a CO transition may oc- 
cur when the inter- site Coulomb repulsion is larger than kinetic energy terms, this is only 
possible in low carrier density semimetals or insulators like a'- NaV205 . Charge order- 
ing can be viewed as a Wigner- crystallization on a lattice®. This should not be confused 
with the CDW transition in more metallic systems. The CO- mechanism in insulating o;'- 
NaV205 can be described within an effective frustraded 2D- Ising model^. It leads to in-hne 
or zig-zag charge order depending on whether the difference in Coulomb repulsion, K1-K2 
between n.n. ( Ki) and n.n.n. (K2) is positive or negative respectively. Later we will also 
discuss alternative CO structures. In Rcf. 8 the possible origin of spin gap formation was 
discussed for the in-line structure where an induced SP transition slightly below the primary 
in-line CO transition was proposed. This scenario would naturally explain the appearance of 
two superposed phase transitions from thermal expansion measurements^ and the observed 
anomalous BCS- ratio. As mentioned the zig-zag CO is an alternative possibility, it has 
been discussed in Ref. 10 and a related structure in Ref. 4. It has been claimed, though not 
discussed in any detail that this structure leads directly to a gap in spin- excitations. 

Important information on the true low temperature CO structure may be obtained from an 
investigation of the complete dispersion of magnetic excitations, especially along a* (_L to 
the chain axis b*). However the existing neutron scattering rcsults^^ were rather limited in 
resolution. A special behaviour of excitations for wave vector q={qx,n){m units of ^ and ^) 
was proposed: The spin gap mode with Ag= 10 meV was suggested to be twofold degenerate 
at q^— 0, 27r and to split into two excitations about 2-3 meV apart for intermediate q^,. This 
was also discussed in a theoretical modeP^. But more recent experiments with much better 
resolution^^ have shown that this is definitely not true and a splitting of ~1 meV exists 
even at the points q= (0, vr) and (27r,7r). Furthermore new electronic structure calculations^"' 
based on the LDA+U approach have shown that there is an additional important exchange 
coupling which has previously been neglected. In addition like in the cuprates small ex- 
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change anisotropies may also lead to gaps for spin excitations. Therefore it is desirable to 
develop a general theory of magnetic excitations in a'- NaV205 that incorporates all these 
aspects and allows to calculate all possible features of the spin excitations in the various 
candidate CO- structures of a'- NaV205 , including the effect of an external field. 

In the following the exchange model for the CO- structures is defined (section II). For the 
low temperature CO phases with intra-chain dimerization it may be mapped to a simplified 
model including only relevant dimer variables (section III). In section IV the spin dynamics of 
various exchange models for a'- NaV205 will be investigated including exchange anisotropies 
and external field. The resulting collective magnetic excitations are studied for all models 
under special emphasis of the importance of intra-chain exchange anisotropies and their in- 
fiuence on the mode dispersions perpendicular to the chain (b-) axis. Finally our calculations 
and their connection to experimental results are summarized in section V. 

II. ELECTRONIC STRUCTURE, CHARGE ORDER AND EXCHANGE MODELS 

In the high temperature phase (T>Tc) a'- NaV205 is an insulating mixed valence com- 
pound whose electronic structure is now reasonably well understood^'^'*. In an effective tight 
binding (TB) model including only V(3d) orbitals one has bonding (B) and antibonding 
(AB) bands corresponding to the symmetric and antisymmetric molecular orbitals of each 
V-V rung. In the following a similar convention for the notation of TB hopping matrix 
elements is used as for the exchange integrals in Figs. 1,6. The B-AB gap ~ i is about one 
eV and the band widths are ~ 0.5 eV (B) and almost zero (AB). This difference has an im- 
portant origin^^ which was not realized previously: Since the dispersion of B and AB bands 
are proportional to t+t^^ and t-t^ respectively it means that tci~ t, and hence t^, the hopping 
across the ladder diagonal cannot be neglected and is necessary for a realistic TB model of 
both B and AB bands. In a naive superexchange model this would also mean that the AF 
exchange constants J~ ^ and ~ ^ should be of the same order of magnitude. This is 
indeed confirmed by spin-polarized LDA-I-U calculations^^ where CO for the 3d- electrons in 
the V-V rungs has to be assumed. They also show that CO a'- NaV205 is in an insulating 
state for sufficiently large on-site U> 3 eV contrary to conventional LDA-calculations which 
predicts a metallic state. As a mean field like theory with broken orbital symmetry the 
LDA-I-U approach does of course not describe the true microscopic nature of the disordered 
MV insulating state above T,.. This is still an open problem. The transition from the high 
temperature MV state to the CO state was investigated in Ref. 8. It was described within 
a frustrated 2D Ising model where the Ising spin r^=±l denotes the 3d- electron localized 
on the right or left position of the rung. In this context the CO of Fig. la and Fig. lb can 
then be described by an order parameter (r^)^ = '^ii'^l) expi^QRi)- For Q—^ one obtains 
the "ferro-" type in- line CO structure and for Q—{7r,7r) the "antiferro-" type zig-zag CO 

— * 

structure depending whether K(0)=Ki-K2>0 or K((5)=K2-Ki>0. In this way model Hamil- 
tonians of the effective Ising type as in Ref. 8 or extended Hubbard models in Hartree-Fock 
approximation^*^ may be used to describe the CO transition in a qualitative way. However 
it is illusory to use such model Hamiltonians in an attempt to actually predict the most 
favorable CO structure. This requires a method hke LDA-I-U which can provide ab-initio 
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(aside from U) total energies of the various CO structures. It has been sucessfuUy used for 
CO- phenomena in semimetalhc 4f- compounds^ and may also be a powerful method for 
the vanadates^"^. In this work the CO mechanism itself is not considered. We rather start 
from plausible candidate structures at low temperature and an appropriate exchange model. 
The derivation of an exchange model for a'- NaV205 from an original extended Hubbard 
model was described in Rcf. 8 and is briefly recapitulated here. It proceeds by eliminating 
high energy charge fluctuations between the rungs, thus confining one d-clcctron or spin 
within each rung. Within this subspace the original model may be mapped to an effective 
low energy Hamiltonian containing the d-electron spin and (Ising) pseudo- spin degrees of 
freedom, the latter describes which of the degenerate V- positions in the rung is occuppied. 
This Hamiltonian is formally similar to those used for the manganites where the pscudospin 
describes an orbital degeneracy of Mn ions. The Ising variable describes the CO transition 
and for T<^ Tc where the intra- rung charge fluctations are also frozen we may replace 
it by its expectation value, i.e. the CO parameter. In this way the effective Hamiltonian 
reduces to an effective spin exchange Hamiltonian only, however with an exchange constant 
Jnm (n,m= V-sitcs) that depends on the CO parameter, i.e. on the degree of charge dis- 
proportionation between the inequivalent V-atoms in a'- NaV205 . In this low temperature 
approximation which we use here the actual size of the CO parameter is absorbed in the 
exchange constants and influences only the energy scale of the spin dynamics, the form of 
the exchange Hamiltonian (T<S Tc) is of the usual type as for spins in a completely CO 
system. In our case due to the orthorhombic symmetry it is essential to include exchange 
anisotropics which may be important for small mode splittings as observed in a'- NaV205 . 
The model exchange Hamiltonian for the proposed CO structures in Figs. 1,6 is then given 
by 

H — - P qxQx I jy qycy , jz oz oz \ 

n,m 

-SI^bHY^S: (1) 

n 

Here J"^ (a =x,y,z) denotes both inter- and intra- chain couplings which may be dif- 
ferent along the three crystal axis a, 6, c (x,y,z). Note that part of the anisotropy in Eq.(l) 
may be due to a Dzyaloshinski-Moriya interaction which can be transformed away in ID in a 
manner described in Ref. 24 and references cited therein. A Zeeman term with field direction 
perpendicular to the Vanadium ab- planes is also included to study the field dependence 
of excitations. Which exchange couplings have to be used depends on the CO structure, 
i.e. the position of the V^^ S=| spins because exchange bonds to V^"*"- ions with no d- 
clcctrons and S=0 are irrelevant for the spin dynamics. This is shown in Figs.l and 6 with 
sets of intra- (J, J^) and inter- chain (J', J'^, J;) exchange parameters (the cartesian index 
a is suppressed). The former may be dimerized to Ji^2=J(l=t^) (Fig. la) and Ji^2=Jd(l±^) 
(Fig. lb). This set has been enlarged as compared to Ref. 8 where only J, J' were included. 
Note that J^^ and J' are not contributing in the in-line CO structure of Fig. la; J is inactive 
for the zig-zag structure (Fig. lb) and J' is not relevant in the structures of Fig.6. In this 
work we consider the following cases: (1) quasi-lD models either in the in-line, zig-zag or 
ladder CO where the inter- chain or -ladder couplings J',J'd etc. are assumed to be much 
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smaller than the intra- chain couplings J,Jd or the intra- ladder J. (2) a quasi- 2D model 
where J' is of the same order as J and J^^. This possibility has been suggested by recent 
LDA+U results. (3) a mixed CO/MV structure which will be discussed later. Different 
methods have to be used for calculating the excitation spectrum in these cases. Figs.la,b 
show that J and Jd play the same role for in-line and zig-zag quasi- ID models respectively. 
Therefore one has in both cases quasi- ID spin chains with intra- chain coupling J (in-line) 
or J(i (zig-zag) coupled by small inter- chain interactions J'^ (in-line) and J' (zig-zag). There 
is one essential difference however: In the in-line structure the CO- transition itself does 
not lead to a dimerization of the chain with an intra-chain J along b. This may be due to 
a secondary SP- transition shghtly below^ leading to a dimerization J — > J(l±5) with S <S 
1. On the other hand CO in the zig-zag structure may itself be accompanied by a lattice 
distortion such that the two legs of the zig- zag chain have different length leading directly 
to a dimcrized exchange Jd{^±6). However it is possible that even in this structure the most 
important contribution to the dimerization comes from the exchange energy along the 
zig-zag legs in Fig. lb. Irrespective of the origin of dimerization a spin gap opens for both 
CO chain structures with its size As(5) depending on dimerization strength. On the other 
hand if a spin ladder structure as Fig. 6a is realised in a'- NaV205 a spin gap will appear 
already without dimerization along b. Finally in the quasi- 2D model with strongly coupled 
chains and in the mixed CO/MV model a broken symmetry spin wave (SW) calculation will 
show that the spin gap can be attributed to pure anisotropy effects. 

For the single dimerized chain or the spin ladder methods based on the Jordan- Wigner 

transformation^'^^'^^ exist to investigate the excitation spectrum. However the focus in this 
work is primarily on the typical behaviour of the transverse dispersion (g ± 6*) of excitations 
where the influence of interchain coupling and exchange anisotropics has to be studied. For 
this purpose it is necessary to use a simple theory as starting point for the intra- chain 

— * 

excitations (|| b*). It is physically appealing to use a spin dimer representation where the 
presence of a spin gap is already manifest in the local dimer basis as a singlet-" triplet" split- 
ting. This representation may also be mapped to the so called bond boson model introduced 
in Ref. 17 for spin ladders. However, for the purpose of investigating spin excitations only 
it is more convenient to keep the original spin- dimer basis, especially when the effect of 
exchange anisotropies on excitations is to be considered. The basic features of the spin- 
dimer representation following Ref. 18 are outlined in the next section and adapted to the 
relevant CO spin structures on the Trellis lattice. 



III. THE LOCAL SPIN DIMER MODEL 

In the dimerized phase of the chain models Fig.la,b or in the case of a ladder with J > 
J (Fig. 6a) it is a useful approach to start from a basis where the strongest exchange pairs 
i.e. J(l+5) or J respectively are diagonalizcd exactly and the weaker couplings are treated 
perturbatively in random phase approximation (RPA). This method has first been used in 
Ref. 18 in a different context. Presently this means the introduction of dimer variables 

Ki — Sii + 82% 
—*—*—* 

Li — Sii — 82% (2) 
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for each pair of strongly coupled dimer spins {Sii,S2i) where Ri denotes the positions 
in the dimer covering lattice. Using this mapping the Hamiltonian in Eq.(l) may be trans- 
formed to 



ia i 

-^EW^J-^ E Js(i^j)LtL'^ (3) 

{ij)a <Cij>a 

Here the first and second term ~ J" describes the local dimer energy and the Zeeman 
energy respectively, the third term ~ J2 denotes the n.n. dimer interactions along the chain 

— * 

direction b and the last term J3 interactions of dimers on different chains. For the two 
chain CO models we have Jf = Ja(l + 5), J 2 = Jq(1 — ^) (with J„ >0 for AF intra- chain 
exchange) and J3 depends on the specific model discussed. Since J and play the same 
role in the in-line and zig-zag model respectively we formally identify J^— > J in subsequent 
discussions of these two models. For the ladder model one has J" —Ja and J"— Jd"J"- 
In the Hamiltonian of Eq.(3) irrelevant parts containing terms ~ K" and L" K" are 
not included because they do not have matrix elements from the singlet ground state to 
the excited states and hence do not contribute to the dispersion of spin excitations^^. The 
energies and states of the S=| dimer are given by 

£;i = 

^2 = Ji = ^(Jf + Jf) 

E, = J^-j[+j, = ^{Ji + Jl) = A' 

li^i) - ^(1 ti) - I it)) (4) 
l^2) = ^(|U) + UT)) 
IV'3) = ^(ITT) + |U)) 
IV'4) = ^{\ tt) - I ii)) 

The ground state singlet {ipi) is separated by an energy ~ Ji from the triplet states 
which are slightly split by an energies j'l-ji and j'l due to the exchange anisotropics given 
by ji=|(J!-Ji)> j'i=i(Ji-Ji) where |ji|, -C |Ji|. In the isotropic case ji= j[ = 
the excited states form a dimer triplet at A = A' = Ji. The dipolar matrix elements 
\M^\'^=\{ipi\La\iJi)\'^ calculated from Eq.(4) are given by |M^p= |M^p=l and zero else. 
Therefore in the dynamical spin susceptibility Uapiuj) of a single dimer only E3 = A' and 
E4 = A appear as possible dimer excitations which is obvious from the form of dimer states 
in Eq.(4). For the present zero field case Ua/3(a;)= Uaa('^)'^a/3 («,/5= x,y) with 
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2A 2A' 

Due to both intra- and inter- chain interactions the two local dimer excitations at A, A' 
will turn into dispersive propagating modes whose minimum energy is the spin gap A^. 
Before we discuss this in detail in the next section we first investigate the effect of an 
external field || c on the dimer states described by the Zeeman term in Eq.(3). Because 
commutes with K'^ one has {tjji\Kz\i'i)=0; i.e. no mixing of singlet and triplet states. The 
only non zero matrix element is {ip^lKzlip^) , therefore the energies Ei 2 and states |'0i,2) will 
be unchanged but E34 and IV'3,4) become field dependent: 

= uliJs) +^1^4) 

= 1 (6) 



The energies of the new eigenstates \%l)±) are given by 

= A'(/i) = J,-j[ + + (1 + {-f)^ 
2 Ji 

E_ = A{h) = J^-j[ + - (1 - (7) 

In the limit h — > |'0±) l'03,4) ^-iid E± — > E3_4. The local dimer susceptibility Ua/3{u;) 
is now given by 

, , 2v?Mh) 2v^A'(h) 



2 



_ 2v?M{h) 2v^A{h) 

. ^ A(/i)2 - A'(/i)2 

UUU^) = -U^M = 2^^^^ (^(;,)2_,2)(A.(;,)2_,2) 

The nondiagonal part is induced by the field. Eq.(8) fully describes the local dimer 
magnetic response and is the basis for the determination of the collective excitations in the 
various CO structures of the Trellis lattice (Fig.l and 6b). 



IV. COLLECTIVE MAGNETIC EXCITATIONS 

In the previous section the effect of the largest intra- dimer exchange interaction J" has 
been treated exactly within the single dimer subspace. The effect of inter- dimer exchange 
may now be treated perturbatively within random phase approximation (RPA). In this 
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method the collective magnetic excitations of the chain or ladder system are given by the 
dynamical RPA susceptibility 

Xiluj) = [1 -Jiq)^u)]-'^Lo) ^D-\q,Lo'MLo) (9) 

Here u(uj) is the local dimcr susceptibility tensor of Eq.(8) and J{q) the exchange tensor 
between the dimers which depends on the specific CO- model considered, q={(lxAy) is a wave 
vector in the reciprocal a*b*- plane in units of ^ and |. The tensors in Eq.(9) have double 
indices: Cartesian a,P — x,y,z as well as CO- sublattice A, r =A,B. Explicitly J^:f(g) = 
(5Q/3J^"(g)and u"^{uj)=SxTUai3{u!)- For two sublattice CO- structures and two local dimer 
excitations A, A' with x,y polarisation one has to expect four (k=1-4) collective excitation 
branches uJk{(1)- They are given as poles of x{q^uj) or zeroes of D(g', c<j). Strictly speaking 
this treatment is only valid when the intra- dimer exchange is appreciably larger than the 
inter- dimer coupling. For example in the dimerized chain models of Fig.la,b the limit 5— > 
is problematic because then Ji — > J2 i.e. intra- and inter- dimer exchange become equal. 
As shown below, Eq.(9) nevertheless leads to the qualitatively correct behaviour for the spin 
gap ^ {)) ^0 although with a different scaling exponent. This indicates that the 

present approach is more effective than the bond- boson theory in MF- approximation^^ 
which leads to a singular A^ for 5= 0. 



A. Excitations for single dimerized chains 

To separate the effects of intra- chain exchange anisotropies from those of inter- chain 
or sublattice coupling it is useful to analyse first the single chain case at zero field. Then 
J^"(g)= J/Tl^l^Ar is diagonal in the dimer sublattice basis and Eq.(9) factorizes for x,y 
polarisation and only two modes exist. The resulting zeroes of Dx^y{q,uj) are then the two 
propagating dimer excitations ujx,y{q) where q = qb* is directed along the chain direction. 
The result applies both for the single linear chain and the zig-zag chain in Fig.la,b (with 
renaming J^; — > J implied as explained in the previous section). Using Eq.(5) and the 
appropriate J{q) the mode dispersions are obtained as 

U;l{q)^l{Jy+J.)[l{Jy+J.){l + Sr 

- J^(l - S'^)cos2q\ (10) 

u;i{q)^l{j.+j.)[l{j.+m+sr 

-Jyil- 6^)cos2q\ 

The spin gap As{S) is obtained as the minimum of u!x,y{q)- The x,y- mode splitting at 
the q=0 is then given by 

ul{Q)-ujl{0) = ^iJ,-Jy)il + 6) 

[(1 + 5)(A + A')-(1-5)J.] (11) 
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which is proportional to the in-plane anisotropy ji= |(J^-Jj,)(l-|-5). If ji = the x,y 
modes are degenerate which can already be seen from their corresponding local dimer exci- 
tations E3, E4 in Eq.(4). For d= one has 

ujx = (AZ),)5 ; = ^[Jy + J, - 2Jx] > 

CUy = {ADy)-^ .Dy^^[jx+J,- 2Jy\ > (l2) 

For the uniaxial case, using >J^=Jy without loss of generality and D2;^j^=|[J2-J2;,j;] 

this leads to ujx,y = \{Jl — Jx,y)^ = ^s- In this limit the spin gap is a pure anisotropy 
gap. Approaching the Heisenberg case Ag vanishes. It is also instructive to consider the 
dispersion Eq.(ll) directly for the Heisenberg case for 5 > 0: 

ujly = uj''{q) = 2J^{l + 5){s\Ti^q + 5cos'q) (13) 

This leads to a spin gap given by 

^ = (2(5)^(1 + (5)^ (14) 

For 5 ^ it vanishes like Ag ~ 5^. Thus the spin dimer RPA approximation gives 
again a qualitatively correct behaviour although the 5- scaling exponent | is smaller than 
the exact one^^ which is |. The dispersion Eq.(13) reduces to 

uj{q) — QjJsin q (15) 

for the undimerized chain whith a — -\/2- This is slightly smaller than the value aocp — 
71/2 for which Eq.(15) describes the lower boundary of the exact Des Cloizeaux Pearson 
(DCP) excitation spectrum of the ID HAF^^. Of course the present spin dimer theory 
completely misses the fact that the excitations really consist of a free two spinon continuum 
since it starts from local dimer excitations which could be interpreted as two spinon bound 
states. 



B. Excitations for weakly coupled dimerized chains, the transverse dispersion 

problem 

The results of the last section give confidence that the basic properties of magnetic ex- 
citations in dimerized spin chains are correctly described by the dimer RPA- theory. The 
advantage of this approach, aside from its simplicity lies in the fact that it can easily be 
extended to include inter- chain coupling. These couplings may lead to transverse disper- 
sion with ql. b*, i.e. a dependence of excitation energy on q^ in addition to the intra- chain 
dispersion or dependence on q^^. As mentioned previously the q^;- dispersion may give im- 
portant clues about the underlying CO- structure. 

First we consider the zero- field case: Then again Eq.(9) factorizes into x,y- polarisations 
but now with sublattice- exchange terms for each polarisation given by {a =x,y) 
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J1a{^ = J%b{^ = Jd{q) 

Jlsio) = Jba{q)* = Jn{q) (16) 

where J^, refer to intra- and inter- sublattice exchange with A,B denoting the two 
inequivalent dimer sublattices of the CO- structures. The four magnetic excitation branches 
of the planar system of chains in Fig.la,b are obtained as solutions of 

l-u^^{u;)[JUq)±\Jm\] = (17) 

The choice of ± in this equation determines the frequency of the acoustical (A) or optical 
(O) mode with respect to the two sublattices. 

In-line chain structure: 
We first discuss the in-line CO structure of Fig. la. It has exchange Fourier components 

Jn{Q) = -JT^^^%^^^\ilx + Qy) (18) 

Together with Eq.(5) the above equations lead to the four mode dispersions 

-Jx{l - S^)cos2qy] 

±{Jy + Jz)Jd\^ + ^) smqysm-{qx + qy) (19) 

^l±{^ = \{J:c + J.)[\{J. + J.)(l + 5f 

-Jy{l - 5^) cos2qy\ 

MJx + Jz)Jd{'^ + smqysm^{qx + qy) 

An interesting aspect of this equation is that for qy= 0, n where the excitation energy 
is close to the spin gap A^, there is no transverse dispersion for the in-line CO model along 
the lines (qa;,0) and (qx,7r). The dispersion of modes for the present case is shown in Fig.2, 
unfolded in the (qxAy)- plane. The inter-chain coupling has its largest effect at the 
maximum mode energy along (q^:,!) where it causes an additional acoustic(A)- optic(O) 
mode splitting connected with the ± in the above equation and in addition it leads to a q^;- 
dispersion. On the other hand when q = (qa;,0) or {qx,'^') J'd has no effect and the observed 
mode splitting in Fig. 3 is dispersionless, it is not of A-0 type but has pure anisotropy 
character as in the single chain case of Eq.(ll). 

zig-zag chain structure: 
In Sec. IV. A it was noted that for a single chain this model is equivalent to the in-line 
structure. However it can be seen from Figs.la,b that the inter- chain coupling is different 
in the two models. For the in-line structure a given dimer is symmetrically coupled with ± 
J'd to four dimers on two neighboring chains whereas in the zig-zag model the coupling is 
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asy metric with strength -J'^, |J'. This leads now to Fourier components for the exchange 
given by 



JD{^ = \j2^os2qy (20) 



13 1 3 1 

Using Eq.(17) we obtain the exphcit solutions (with the formal replacement —>■ J) 

liJy + Jz)[liJy + Jz){l + - 7.(1 - 5') cos2qy] 

3 11 3 1 

±[{Jy + J,)[J'f cos{-qy - -q^) - -J'^'cosi-qy + -q^)] 

^1±(q) = (21) 

liJx + Jz)[liJx + Jz)il + Jy{l - 6') cos2g,] 

3 11 3 1 

±[{Jx + Jz)[JdCOs{-qy - -q^) - -J'^cos{-qy + -q^)] 

While the intra-chain part in this expression is the same as in the in-line model of Eq.(ll) 
the second part leading to the transverse q^;- dispersion is completely different. For example 
taking qy— tt we obtain 

u;l_{q^,7r) -a;^+(g^,7r) = 
2{J'i + ^J'n{Jy + Jz)sm^q. 

ujl_(qx, tt) - u;y_^_(q:c, tt) = (22) 
2{f,' + ^J'niJx + Jz)sm^q. 

This shows that in addition to the anisotropy induced x,y- mode splitting each of them 
shows a further splitting into A,0 (±)- modes which has dispersion: it vanishes at qx=0 
and is at maximum for qx='n'. This situation is clearly illustrated in Fig.3. Whether this 
dispersion is visible in the experiment depends on how large it is against the pure anisotropy 
splitting caused by the intra- chain exchange. In principle both are present and Fig. 3b shows 
two typical possibilities. The qa; dispersive A-0 splitting which is absent for qy= 0, vr in the 
in-line case therefore in principle offers a possibility to distinguish between both models. 

Finally we discuss the intensity variation of qy= 0, vr spin gap modes as function of total 
momentum transfer k= q + f [q & l.BZ) mentioned in Ref. 11. It was observed that the 
intensity of the hw— 10 meV excitation exhibits unexpected variation in with period h=3 
where t— (27rh,27rk,0) is a reciprocal lattice vector in the ab- plane. For a strictly ID system 
the intensity should rather be constant and therefore this variaton possibly points to a more 
2D character of magnetic excitations. We now analyze the intensities in the dimer RPA 
model for that structure. For simplicity we neglect the additional splitting of modes caused 
by xy- exchange anisotropy, i.e. we assume ^x=^y Then the observed splitting along q^; is 
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entirely an A-0 splitting due to the fact that the dinner lattice consists of two sublattices. In 
this case the intensities may be obtained from the dynamical susceptibilities^^ decomposed 
according to 

X(g + f,a;) = ^(l + cos$)x^(g,o;)+ (23) 
^(1 -cos$)xo(g,^) 

where (j)=f ■ p and p={\-i -\) is the vector joining the two dimer sublattices in units of 
a,b respectively (Fig. lb) which leads to $=h7r-k7r. Prom the imaginary part of x(g, a;) the 
A,0 intensities are obtained as {q e l.BZ) 

Ia,o{Q + ^) = 7r^T^[l ± cos(/i7r - kTT)] (24) 

where ± corresponds to A,0 respectively. In the experiments^^ one has q={(lx,'^) and r 
given by (27rh,0). Neglecting the small A-0 splitting, i.e. setting u!a,o{q) — one then has 

/A,o(r) = ^(l±cos/i7r) (25) 

Two points are worth noting: The period of the intensity is given by h=2 and not h=3. 

The intensity maxima of slightly split A,0- modes are shifted by one half period (h=l). 
Experimentally the intensity at an energy transfer hw= 10 meV was measured as function 
of h. Since this energy is just in between upper (O) and lower (A) mode and both have a line 
width considerably higher than their splitting the measured intensity is then the average of A 
and O mode intensity. According to Eq.(25) however the average is a constant independent 
of h, irrespective of the period of individual A,0 intensities. We conclude that the zig-zag 
CO structure, at least in the dimer RPA model for weakly coupled zig-zag chains, does 
neither explain the observed intensity variation nor its period. 

field dependence of excitations: 
Investigation of the field dependence of magnetic excitations may give further information 
on the nature of the spin gap and its observed splitting and transverse dispersion. The 
field dependence may also be calculated from the basic Eq.(9) where it enters through the 
local dimer susceptibility Eq.(8). Due to breaking of time reversal symmetry there is now a 
nondiagonal term Uxy= -u*^. and both polarisations couple to each of the A(h), A'(h) local 
dimer transitions. The poles of Eq. (9) are then given by 



1 - [Uxx{i^)J±{q) + Uyy{uj)jl{ 

+ [uxx{uj)uyy{uj) - \uxy{uj)\VUQ)J±iQ) = (26) 

Here J"(^= J^(^±J^(^ and ± has to be taken synchronously at all positions. After 
straightforward but lengthy algebra the solution of this equation leads to the field dependent 
dispersions for the magnetic excitation branches uj^iq, h) {k =1-4): 



12 



C±{q,h) 



^B^±^[Bl{q,h)-4C,{q,h)f^ 
A\h)+/\'\h) 

-2[A{h)Jl{q) + ^'{h)Jl{q)y 
A\h)A'\h) - 

-2[A{h)Jl{q) + A'{h)Jl{q)]A{h)A'{hy 
-2[A{h)J^{q) + A\h)JUq}]Aih)A'{h)v^ 
+4J^{q)Jl{q)[{^\h) + A'2(/,))mV + 
A(/i)A'(/i)(M^ + t;^)] 



(27) 



Here a = ± and k, = (±, a) corresponds to any of the four possible combinations of ±- 
signs in the last equation. In the quantities B± and C± the ± signs always have to be taken 
simultaneously. With u(h) and v(h) given by Eq.(7) and A{h), A'{h) by Eq.(7) the above 
expressions represent the complete solution for the field dependent dispersion of magnetic 
excitations in the anisotropic coupled dimcr system. These equations can be applied to the 
CO structures of Figs. 1,6a. The specific CO determines only the exchange functions 3%{q)- 
For zero field this equation reduces to the previously studied solutions of Eq.(17). Figure 
4 shows the field dependence of g=(0,7r) modes, i.e. the spin gap modes vs. external field 
for the zig-zag CO in the two limiting cases corresponding to Fig. 3b. One obtains a quasi- 
linear Zeeman splitting of g=(0,7r) modes in the small anisotropy (2) case and almost field 
independent modes for large anisotropy (1). The gap will close only at a very high field 
which is expected since the zero field spin gap of A^ =10 meV is quite large. In this model 
it was assumed that the dimerization S itself shows httle field dependence since it should be 
a direct consequence of the lattice superstructure induced by the CO. 
As in the zero field case the susceptibility a;) in Eq.(9) may also be used to calculate 
the intensity of the four uiK{q, h) excitation branches. They are given by the imaginary part 
7r~^X^(g, a;)" where X^=\{Xxx + Xyy)- obtains delta- function contributions of the type 
Zi^5{uj — uji^^q)). The intensity of each mode u)^{q, h) can be obtained from Eq.(26) as 



In the isotropic case with J^= Jj_ and A{h) = A'{h) this reduces to a simple formula for 
the two (A,0) modes {a — ±): Z^{q,h.)— A{h)/ijj^{q, h) which corresponds to the prefactor 
of the zero field result in Eq.(25), the variaton with r is supressed here. Within the field 
range of Fig.4 there is only a few per cent change of the corresponding mode intensity. 



One reason for focusing on ID chain models for the magnetic excitations of a'- 
NaV205 was the observation of quasi- ID temperature dependence of the susceptibility 




(28) 



C. Excitations in strongly coupled chains 
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in the MV phase above T^.. This was attributed to d-electrons locahsed in the molecular 
bonding orbitals of each V-V rung having strong exchange J along the ladder and weak ex- 
change J' between them. This picture was qualitatively supported by by LDA- calculations^ 
mapped on an effective 3d- tight binding (TB) model which lead to very small hopping 
matrix elements t'<St suggesting that <^ '^~ir ^ simple superexchange picture. 

However a recent LDA+U analysis"*^^ with a mapping to an extended TB- model including 
both V3d and 02p orbitals has seriously questioned this picture for the low temperature 
CO phases. In this calculation the mapping of LDA-I-U total energies of various CO and 
spin polarized states to that of a corresponding Heisenberg model enables one to calculate 
realistic values for the most important exchange constants. It turns out that in the CO 
phase J' is only about a factor of two smaller than and this "diagonal" ladder exchange 
is even bigger than the exchange J along the leg of the ladder. Furthermore surprisingly 
even the J'^ exchange constant is not negligible and both J' and J'd are ferromagnetic. For 
a realistic value of U= 3eV the exchange constants have values as given in the caption of 
Fig. (5). The reason for the large J' in the CO structure as compared to the homogeneous 
MV state lies in the change of pd- hybridisation due to the shift of 3d- levels on inequivalent 
V- atoms^^. If this LDA-I-U result for the exchange corresponds to the real situation then 
CO a'- NaV205 is magnetically more like a 2D system with strong AF coupling along the 
ladder diagonals and legs and almost equally strong FM coupling between the ladders. Such 
a model is very different in principle from the ID models discussed sofar. We now also 
investigate its magnetic excitations and origin of the spin gap which is different from the 
dimerization mechanism in this model. This is also partly motivated by the fact that accord- 
ing to Ref. 21 there is indeed no intra-chain dimerisation in the low temperature structure 
as assumed in the previous models. Naturally the dimer approach of previous sections is 
not possible for the zig-zag CO structure with its very large interchain coupling J'. On the 
other hand there is no problem for the in-line structure since J' is inactive in this case and 
even the appreciable J'^ obtained from LDA-I-U does not affect the dispersion very much 
since it is effective only at the maximum energy and docs not influence the spin gap as 
shown in Fig.2. For the zig-zag CO instead we now start from a broken symmetry ground 
state with a spin configuration as indicated in Fig. (lb) which has the lowest ground state 

T' J' T 

energy E= y — ^ — 4 ■ Of course this approach does not describe the real ground state of 
a'- NaV205 which is nonmagnetic, nevertheless the excitation spectrum can be expected to 
have realistic features. The spin state consists of four magnetic sublattices Af , Bf , AJ,, BJ,. 
The molecular field for the sublattices A is given by A(|)=-A(J,) = A with 

A^{S)[J''-2J'i-2r,] (29) 

As in the previous models we include anisotropics in the largest exchange although 
its magnitude has not yet been calculated in LDA-I-U which was applied without spin- orbit 
coupling^''. Without loss of generality we assume that the spins are oriented along the c-axis, 
i.e. >Jd'^- Furthermore (S) is the saturation moment equal to 1/2 at T=0. The RPA 
equation for the spin wave (SW) modes is formally the same as Eq.(9) but the dynamical 
variables are now the individual spins and not the dimer excitations. Therefore instead of 
Eq.(5) for the dynamical suscepibilities we have now Uxx{i-^)= '^yy{^)= ^(0;) and Uya;(a;) = 
Ua;j/((^)* = v(a;) with 
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where {S)x—±{S) and Aa=±A for A =t) i sublattices. Furthermore the exchange 
Fourier transforms JoNio) i'^ =x,y,z) are now tensors defined in the original spin lattice 
(1,1 sublattices) instead of the dimer covering lattice. The various components connecting 
the four sublattices can be read off from Fig. (lb), e.g. Jd|| =J'(^) ^tc. with 

J'{q) = J'expi^i^q^ - Qy) 

J^(g) = 2J^cos^(g, + g,) (31) 

Jd{q) = 2JdCxp{-^q^)[cosqy - i6 sin qy] 

Here 5 ^ 1 is the dimerization of the zig-zag chain which may exist due to the low 
symmetry of the corresponding CO structrue. After some algebra the complete RPA spin 
wave solution of Eq.(9) applied to the present case consists of four branches (k =1-4) which 
are given by cUk— ^k/ (S) with 

±{{c[±d[)^ - [|ci|2 + - |C2|2 - \d2\^ 

±{cidl + cldi - C2dl - cld2)]}^ (32) 

ci = A^ + j'^j'y* - {j'fj'y + r.jy;) 
C2 = -A(j- + j'y) - {r.j'y + jpf) 
d, = A{j'f - jf) + rjf - j'y*Ji 
d2 = A{Jl - Jl) + rj'! - f'J'i 

Here A=A/(5') and c'i=(ci-Fc^)/2, d'i=(di-Fd^)/2 denotes the real part of these func- 
tions. Note that the variable q was suppressed in J'(^i J'd(^ and Jd{q) in the above 
expressions for simplicity. The signs ± with a hat have to be taken simultaneously with 
upper or lower value wherever they appear thus leading to four spin wave branches. Again 
it is useful to consider the solutions for the single chain case only, i.e. setting J'=J'ci=0. 
Then Eq.(32) reduces to u;l—c\ ±1^2! or 

coliq) = [A±2.r,^i [At2 Jhd (33) 
7q- = (cos^ qy + (5^ sin^ qy) 2 

Each branch is twofold degenerate since there is no A-0 splitting without inter- chain 
coupling. At zero wave vector one has for >0, using 2(5') =1: 

u;+(0) = [Jj - JS]Hr, + J|]i {2fjaf^ 

0040) = [Jj + J,1^[JJ - J|]^ {2jp,f^ (34) 

with j^'^= J^- J^'^ and Jd=(Ji^+J^+J^)/3 denoting the exchange anisotropy and average 
respectively. 
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This shows that in the SW approximation the spin gap A"=(2jJ| J^;)^ is entirely an 
anisotropy gap and independent of CO induced dimerization. Indeed in the isotropic case 
uj{q)=2{S) Jd(l — 5^)2sinq'j/. The dimerization does not remove the gapless excitations but 
only changes slightly the spinwave velocity. Therefore not surprisingly for isolated dimerized 
isotropic HAF chains the SW RPA approximation is quahtatively incorrect and the dimer 
RPA approach of previous sections should be used. However when interchain couplings J', 
J'rf become larger than J^^ the situation is reversed and the SW approximation of Eq.(32) 
is a better starting point for the essentially 2D magnetic system. This is certainly the case 
when one uses the exchange parameters obtained from LDA+U where due to \J'/Ja\ — 0.5 
J' is much larger than JS with a 5=0.034 estimated from the spin gap in the dimer model. 
The SW excitation branches for the 2D exchange model as obtained from the LDA+U 
parameters are shown in Fig. (5). The size of the anisotropy which is not determined by 
LDA+U is fixed by the size of the spin gap of ^ 10 meV as suggested by Eq.(34). For zero 
anisotropy one would get a Goldstone mode (A- branch) also for the 2D model. Thus the 
anisotropy and not the CO induced dimerization of the exchange is the origin of the spin 
gap in this 2D model. Fig. (5) exhibits a pair of two strongly split A,0 modes where the 
splitting is approximately given by [J' {J' — 2Jci)]2. Both modes show a small additional 
splitting caused by the xy- exchange anisotropy of J^; if it vanishes the A and O modes 
are twofold degenerate throughout the BZ. This fact is well known already from the simple 
(two magnetic sublattice) AF where only A modes exist and the degeneracy can simply 
be understood as a result of the downfolding into the AF BZ. The A- modes are nearly 
dispersionlcss along q^^ because the effect of the AF and the FM J'^ nearly cancel along 
this direction. The behaviour of the A-modes around the (0,7r)- point with their splitting 
increasing towards (7r,7r) is qualitatively very similar to the two modes in the dimer model 
with weakly coupled chains (Figs. 2,3). Note however that the role of A-0 sphtting and 
anisotropy splitting are reversed. It is the latter which now leads to the dispersion of the 
small A-mode splitting whereas the A-0 splitting is now much larger. The existence of the 
high energy split- off O- branch is an essential prediction of this 2D model and could be 
tested directly experimentally. The ffat part of the O- branch lies at about an energy given 
by ujo — [J' {J' — 2 Jd)] 2 ~ 30 meV roughly twice the maximum energy investigated so far^-*^. 

D. Excitations in Ladder and partly mixed valent structures 

Recently a model for the lattice distortion caused by the charge ordering below Tc has 
been proposed based on low temperatue x-ray scattering^^. According to this model only 
every second ladder in Fig.lc is distorted perpendicular to the chain (6-) direction with a 
period of 2a along the (a-) direction. It is not immediately obvious which CO structure 
is compatible with this distortion pattern but the simple in-line and zig-zag model cannot 
easily be reconciled with the curious fact that every second V-ladder of the Trellis lattice 
is undistorted. Two possible models which incorporate this fact have been proposed and 
are shown in Fig. 6: (i) CO structure of the spin ladder type^^'^^ where two V(3d) electrons 
occupy every second rung and the other rungs have no V(3d) electrons. This ladder model 
is therefore completely different from the chain models which have only one 3d- electron 
for every rung. It seems rather surprising that the ladder structure should be realized since 
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LDA+U calculations^'^ indicate that it has a much higher total energy than the chain struc- 
tures, (ii) structure with alternating CO zig-zag chains and disordered chains^'^. Here half 
of the ladders have zig-zag CO like in Fig. lb but the other half remains in the disordered 
mixed valent state. This is again a structure with only one d- electron per rung on the 
average and therefore it will have a total energy not too different from the CO structure of 
Fig.lb. 

We discuss first the magnetic excitations in the ladder structure in qualitative terms assum- 
ing an AF superexchange J" of spins within the ladder rungs via the intervening oxygen. In 
this structure there is no connection between the spin gap and the doubling of periodicity 
along b since it appears already in the undimerized (equidistant) ladder. Essentially one 
deals with single ladder excitations in this case because the magnetic S=| V^"*" ladders are 
separated by nonmagnetic S=0 (V^"^) ladders. There is another important difference to the 
chain structures: The doubling of the periodicity of the lattice (period 2b) does not show 
up in the spin ladder which still has periodicity b as seen in Fig.lc, this has drastic conse- 
quences for the excitations. In the present case the local dimer susceptibility is determined 
by the states of the rung- dimer with interaction J" (Fig. Ic) and it may be obtained by 
replacing A j{J^ + J^) and A' j{J^ + J^) in Eq.(5). Using Eq.(17) which also holds 
for the present case we obtain for the two excitation branches (there is no A-0 splitting in 
this case): 

^lio) = \{Jx + Jz)[\{Jx + Jz) - 2J|cosg,] (35) 

The most important aspect of these excitations is the doubling of the periodicity 27r 
along the ladder instead of tt for the chain models. In the latter the points qjy=0 and ^y—'K 
are degenerate and their excitation energy is at (or close) to the minimum, i.e. equal to 
the spin gap (Figs. 2, 3). On the other hand for the ladder structure assuming Jg=J2-J" < 
the excitation energy is equal to the spin gap only for ({y='K., whereas it is at the maximum 
in the zone center (qj/=0) for AF effective inter-(rung) dimer exchange Jg. The possibility 
of the spin ladder model can therefore in principle be directly experimentally investigated. 
Finally the splitting along the transverse q^^ direction will be constant and determined by 
the exchange anisotropy as in the in-line chain case. 

The more recent proposaP^ of an alternating CO/MV low temperature structure of a'- 
NaV205 which is shown in Fig. 6b is a very interesting possibility and deserves a detailed 
analysis of its magnetic excitation spectrum. In this structure CO zig-zag chains alternate 
with ladders in the MV state along the a- axis. In both we have one d-electron per rung 
on the average but in the MV ladders the electron is not localised on one of the two V- 
positions of the rung but resonates, i.e. it is in the molecular bonding state of the two 
V-atoms. This means there are three types of V-sites with formal valencies Z=4 or 5 on the 
CO zig-zag chain and Z=4.5 on the MV ladder. Whether the existence of three inequivalent 
V- sites is compatible with NMR results is not clear. To describe the magnetic excitations 
in such an inhomogeneous state we make a drastic but reasonable assumption: Since the 
d-electrons residing in the molecular orbitals are spread out over the whole rung, their spin 
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response will be concentrated around zero total momentum transfer contrary to the atomic 
spins on the CO zig-zag chains where the magnetic scattering intensity varies with the 
atomic form factor. Therefore the contribution of the spins in the molecular orbitals to the 
scattering cross section should be negligible at large momentum transfer and consequently 
only the atomic spins residing on the zig-zag chains will be included in the model exchange 
Hamiltonian for the strucurc of Fig. 6b. This model Hamiltonian has only two (in general 
anisotropic) exchange constants: intra-chain similar as in the zig-zag model of Fig. lb and 
inter-chain coupling Jf which connects V"^"*" spins on next nearest ladders via a superexchange 
path across the intervening MV (non- magnetic) ladder. Note that in principle the lattice 
distortion in which was described in the beginning of this section implies a dimerization 
Ji,2= Jd(l=t^) of the intra- chain along the a-direction. This means that the intra-chain 
exchange within a given zig-zag chain is Jrf(l+5) and 3 4(1-6) on two adjacent zig-zag chains 
separated by a MV ladder (Fig. 6b). This type of dimerization therefore leads to a doubling 
of the unit cell of the exchange Hamiltonian along a whose consequences we will discuss 
later. Note that it does not lead to an alternating exchange within a given zig-zag chain 
and hence not to a dimerization gap in the excitations of the isolated chain. For this reason 
it is possible to start from a Neel ground state and use the spin wave approximation even 
though we consider now weakly coupled zig-zag chains. Formally the SW- calculation is very 
similar to that performed in section IV. C and we do not repeat the details. The Fourier 
components of the exchange entering the dynamical susceptibility x(g, cu) may be read off 
from the structure in Fig. 6b and are given by 

1 

Jdio) = cosgyexp(-g^) 

2i 

J^{q) = 2J^ cos Qy exp{—qx) (36) 

The spin wave branches are then again obtained from the poles of x(g, uj) which leads to 
the secular equation detx(g, a')=0. In this equation the dimerization d{\\ a) appears only in 
0(5"^) and if these terms are neglected for the moment we obtain the four spin wave branches 
in the small BZ (| q^ | < f ) as 

= [(Jf - JdJ!) + m sin^ %\ ± Jl^Jl - J!) cos qy 

±Ji{J2 + Jl) cos^ qy cos qx (37) 

Here ± signs are chosen independently and therefore k =1-4. For the isolated chain with 
uniaxial exchange (J/=0, J^^J^ = Jd) this reduces to a single dispersion 

u;'(g^) = + J^sin^g^ (38) 

with the Ising anisotropy spin gap A^= (J^^-J^) — 2Jd(Jd-Jd)- In the Heisenberg case 
Aa=0 and one obtains the dispersion Eq.(15) but now with asw—^ which is smaller than 
the values of a in the dimer approximation and in the exact result (see below Eq.(15)). For 

Jrf 7^ "^l there are two anisotropy gaps A^. 

We now focus on the most interesting part of the transverse dispersion along q^; of the 
coupled chain system. For qy=0 (or qj;=7r, this leads only to an interchange of modes) one 
obtains in the large BZ ( | q^ | < tt) two unfolded modes given by 
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— S'^ cos q. 
A~ — S" cos q. 



'X 



\x 



(39) 



with A±2 = 2JJ(JJ - JY) and 5± - {JdJi)/^'^. In deriving Eq.(39) from to Eq.(37) we 



assumed that 5^ <Sl. These spin wave dispersions are equivalent to the recently proposed 
empirical dispersion relations^^ obtained by new inelastic neutron scattering results which 
had much higher resolution than those performed in Ref. 11. The above formulas provide 
an excellent fit to the experimental dispersions as shown by Rcgnault et al in Ref. 22 with 
parameters obtained there as A+= 10.65 meV, A~= 8.75 meV, 5+= 0.4 meV, 5^= 0.5 
meV. Note that this dispersion proposed by Regnault et aP^ and derived in the present 
model has only half the period in q^. as compared to the pure zig-zag model of Fig. lb 
and Eq.(22). Using the empirical parameters from above and Eq.(39) one may completely 
determine the relevant parameters in the present theoretical model: J^= 38 mcV, Ji= 0.11 
meV for intra- and inter- chain exchange respectively and for the intra- chain anisotropics we 
have ^z-^x— 1-49 meV and ^z-^y— 1 meV. Note that our model result for 5'^ below Eq.(39) 
requires /8~= A+/A^. Experimentally /5~= 0.80 and A+/A~= 0.82 which is close to 
expectation. The dispersion obtained from Eq.(37) (or approximately from Eq.(39)) with 
these parameters is shown in Fig. 7. Finally we comment on the question of intensities. 
Since we have two chemical sublattices separated by li = |a+6 (Fig.6b). and only sublattices 
of opposite spins are coupled in pairs the intensities will be given by a similar expression as 



In this case then the period of the intensity variation is correctly given by h=3 as observed 
experimentally contrary to the zig-zag model of Fig. lb. 

Sofar the dimerisation 5 along a leading to — > (1±5) J^; has been neglected. If included it 
would lead to an opening of an additional gap at the points qa;=±| of size 2(5A(j as shown 
in Fig. 7. Sofar this dimerization gap has not been identified. It is not clear whether the 
resolution is still too small or whether the inter-chain exchange dimerization J(i(l±5) of next 
neighbor zig-zag chains as shown in Fig. 6b is indeed negligible. 



We have proposed and analyzed a number of spin-excitation models that may be rel- 
evant for inelastic neutron scattering investigations in the CO low temperature phase of 
a'- NaV205 . The most frequently discussed models are based on in-line or zig-zag chain 
CO and assume that the exchange coupling along the chains is much smaller than between 
the chains. In these models there is a strong dispersion along the chain axis h* caused by 
the exchange coupling along the legs (in-line J) or ladder diagonals (zig-zag J^). Recent 
LDA-I-U calculations^^ suggested that both are of the same order of magnitude. Therefore 
even in the completely CO zig-zag structure there is a strong dispersion || h* {%)■ The 



in Eq.(24): 



1 Itt tt 

lA{T)^-{l + COs{-h))^COs\-h) 

/o(f) ~i(l-cos(|^/.)) = sin^(|/i) 



(40) 



V. SUMMARY AND CONCLUSION 
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minimum excitation energy, i.e. the spin gap in this scenario is due to a dimerization 
J(l±5) or J(i(l±(5) of the intra- chain exchange. The dispersion || b* (q-^) on the other hand 
is comparatively smalL However it shows characteristic differences for the two quasi- ID 
models. Allowing for small anisotropies in the largest exchange J or we find that (1) the 
in- line model has a q^;- dispersionless splitting of the spin gap mode for qj,= 0,7r determined 
by the exchange anisotropy alone. The inter- chain coupling may only contribute to the q^;- 
dispersion at the maximum mode energy at qj/=f • This is a direct consequence of the Trellis 
lattice structure with every second ladder shifted by |. (2) in the zig-zag model the splitting 
of the qj,= 0,7r spin gap modes has both a contribution from exchange anisotropy J^-J^ and 
inter- chain coupling J'. Depending whether the former or latter is stronger one has little 
or noticeable dispersion along q^ and the role of anisotropy splitting and A-0 splitting are 
interchanged. It was also shown that the magnetic field behaviour in the two limiting cases 
of the zig-zag model (Fig.4) is different. While for J' appreciably larger than Jjj-Jj^ there 
is an almost linear Zeeman splitting of spin gap modes, they are almost field independent 
for small fields in the opposite case. In addition we discussed the zig-zag CO model in the 
case of strong coupling between the zig-zag chains because LDA+U calculations predict a 
surprisingly large inter- chain coupling J' in this CO- structure. In this case we used a 
broken symmetry approach to calculate the spin excitations. In this model the spin gap is 
a pure anisotropy gap, furthermore one obtains a split- off optical branch at an energy of 
30 meV. The observation of such a mode would be crucial for this model, sofar there is no 
experimental evidence that it exists. 

Furthermore we briefly discussed the alternative spin ladder model of Fig. 6a with AF cou- 
pling in the ladder rungs. It was found that the dispersion along q^ has twice the periodicity 
as compared to the in-line and zig-zag chain models. 

Perhaps the most promising model investigated is the mixed CO/MV model with zig-zag 
chains separated by disoredered (MV) ladders. This structure model has been proposed in 
Ref. 23 and we have shown that it leads to spin wave dispersions exactly as those empir- 
ically proposed by Rcgnault et al^^ from new inelastic neutrons scattering results. Most 
importantly it shows half the period for the dispersion along q^. as compared to the other 
models and also has the proper intensity variation. In this model the spin gap is due to 
a predominately Ising type exchange anisotropy and the dimerization perpendicular to the 
chains has only little effect. 

All models discussed account for the basic qualitative properties of the available neutron 
scattering results^^: (1) Sligthly split spin gap modes at As= 10 meV with little or no 
dispersion along q^. and (2) a large dispersion of magnetic excitations along the chain (q^)- 
direction. What is different in the models is the interpretation of the origin of various gaps 
and splittings observed, i.e. whether they are due to dimerization, anisotropy, ladder type 
or of A-0 nature. To make further progress in discriminating between these models (and 
possibly others not investigated here) and also to obtain a more reliable set of exchange 
parameters for them it is necassary to have inelastic neutron scattering results in a larger 
energy and momentum region and with enhanced resolution and also a more detailed in- 
formation on the momentum dependence of the intensity for each individual mode. The 
investigation in this paper has given a clear classification of the typical signatures of the 
different exchange models one has to look for. 
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FIG. 1. Charge ordered (CO) structures of a'- NaV205 discussed in the text. Hatched circles: 
(8=5), open circles: V^"*" (S=0). Oxygen atoms on the legs and rungs of V- ladders are not 
shown. Thick lines (Ji or J) denote the dimer basis of each model, (a) in-line CO with active 
exchange constants Ji,2= J(l±(5) and J'^^; 5= dimerization strength along b. (b) zig-zag CO with 
active exchange constants Ji,2= Jd(l±(^), J' and J'^^. The ib signs denote the spin configuration 
with lowest energy for the LDA+U exchange parameters. This is only relevant for the 2D spin 
wave scenario of Sec.IV.C 
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FIG. 2. Spin excitations 0,i^{q) = u^ictj/J^ (k=1-4) in the in-line structure calculated with 
dimer RPA- theory of Eq.(19). Only J" (a=x,y,z) and J'^^ are active exchange constants, the 
former determines the large dispersion along h* (q^), the latter the A-0 splitting Ayio- The spin 
gap A5 is mainly caused by the dimerization 5. J'^ has no influence along paths with qj,= 0, tt, 
therefore in this model is a pure anisotropy splitting determined by Ja;-Jy and the two split 
modes are dispersionless along q^;. In the xy- isotropic case A^ =0. Exchange parameters used are 
J^= 38.4 meV, Jj,= 37.4 meV, 3^= 37.9 meV, J'd= -6 meV and dimerisation 5= 0.034 
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FIG. 3. Spin excitations 0.^{q) = ujfi{q)/J^ (k=1-4) in the zig-zag CO structure according to 
Eq.(19) (where the formal replacement — > J was made). Only J^, J' and J'^ (Fig.lb) are active 
exchange constants in this structure. (a)In addition to the anisotropy splitting along q^^ there 
is a further A-0 splitting superposed. Which one is more pronounced depends on the relative 



size of xy- anisotropy J^-J^ and the inter- ladder coupling J'. In this plot J^- 



=38.2 mcV, J^=37.6 



meV, J^=37.9 meV, J'= 0.5 meV, J'cZ=0 and 6= 0.034 was used, (b) Enlarged excitation branches 
between (tt, tt) and (0,7r) for two extreme cases, above: large anisotropy JJ^=38.2 meV, J^=37.6 
meV, J^=37.9 meV and small J'= 0.1 meV. below: small anisotropy J2=37.95 meV, J^=37.85 
meV, J^=37.9 meV and large J'=0.5 meV. Other parameters in both cases as in (a). It is seen that 
anisotropy splitting and A-0 splitting A^q interchange roles in the two cases. The situation 
proposed previously^^ corresponds more to the lower part. Parameters for the dispersion in (a) are 
between these two extreme cases. 
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FIG. 4. Field dependence of modes at 
^=(0,7r) for the zig-zag CO. Exchange param- 
eters for (1) and (2) are identical to those of the 
upper and lower part of Fig.Sb respectively. 
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FIG. 5. Spin excitations for the 2D exchange 
scenario described in Section V. The SW results 
of Eq.(33) has been used. The exchange param- 
eters are basically those from LDA+U calcula- 
tions (U=3 eV)^^ for the spin polarized zig-zag 
CO structure of Fig. lb, except for a slightly 
larger J',^ and the additional anisotropics which 
have been introduced to obtain spin gap 
and anisotropy splitting Aa- Explicitly wc use 
j;^=34.5 meV, J^=35.2 meV, J^=36.2 meV, J'= 
-17.8 meV, J'rf= -6 mcV. In the SW picture the 
dimerization does not lead to a spin gap and 
was set to zero. is mainly determined by the 
(xy)-z anisotropy and A^ by the in- plane xy- 
anisotropy. The large Aaq- gap is caused by 
the large J' in this parameter set. 
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FIG. 6. (partly) charge ordered (CO) struc- 
tures of a'- NaV205 discussed in Sec. IV. D. 
Hatched circles: (S=^), open circles: V^"*" 

(S=0). Hatched ellipse: MV state of V-V 

rung. Oxygen atoms on the legs and rungs 
of V- ladders are not shown, (a) ladder CO 
consisting of isolated S=^ ladders with active 
exchange constants J, J and J^. (b) partly 
zig-zag CO/MV structure. Intra-chain exchange 
Ji,2= Jd(l±5)) is dimerized perpendicular to the 
chain under the assumption of the low tempera- 
ture distortion pattern (arrows) proposed in Ref. 
21. Next neighbor chains are coupled by J;. 
The hatched ellipse denote one d-electron in the 
molecular state of the rung, i.e. the V- atoms 
have formal valence 4.5. 
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FIG. 7. Spin wave dispersion (^^(g) along 
q^; for qj^=0,7r unfolded in the large BZ corre- 
sponding to the undistorted structure. Gap- 
less curves: Modes for zero dimerisation ac- 
cording to Eq.(37), gapped curves: modes with 
perpendicular dimerization according to Fig. 6b 
with 5=0.01. This leads to small gaps at the 
boundary (±|) of the small BZ corresponding 
to the distorted structure. Exchange parameters 
and anisotropies have been choosen to comply 
with the experimentally determined dispersion 
by Regnault et al [22] and given in the text of 
Sec.IV.D. 
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